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We study three dimensional quenched Quantum Electrodynamics in the bare vertex approxima- 
tion. We investigate the gauge dependence of the dynamically generated Euclidean mass of the 
fermion and the chiral condensate for a wide range of values of the covariant gauge parameter £. 
We find that (i) away from £ = 0, gauge dependence of the said quantities is considerably reduced 
without resorting to sophisticated vertex ansatze, (ii) wavefunction renormalization plays an im- 
portant role in restoring gauge invariance and (iii) the Ward-Green- Takahashi identity seems to 
(N ' increase the gauge dependence when used in conjunction with some simplifying assumptions. In the 

Landau gauge, we also verify that our results are in agreement with those based upon dimensional 
rcgularization scheme within the numerical accuracy available. 
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I. INTRODUCTION 



In a gauge field theory, Green functions transform in a specific manner under a variation of gauge. These trans- 
formations carry the name Landau-Khalatnikov-Fradkin (LKF) transformations, These were also derived by 
Johnson and Zumino through functional methods, |^|,[|. As a consequence of gauge covariance, Green functions obey 
certain identities which relate one function to the other. These relations have been named Ward-Green- Takahashi 
identities (WGTI), l^). At the level of physical observables, gauge symmetry reflects as the fact that they be 
independent of the gauge parameter. Perturbation theory respects these requirements at every level of approxima- 
tion. However, this has not been achieved in general in the non-perturbative study of gauge field theories through 
Schwinger-Dyson equations (SDE) carried out so far although significant progress has been made. The gauge tech- 
nique of Salam, Delbourgo and later collaborators, [p|-|l4||, was developed to incorporate the constraint imposed by 
Q_i| WGTI. However, as pointed out in Jl^] , gauge technique can become completely reliable only after incorporating 
JL . transverse Green functions with correct analytic and gauge-covariance properties. Another method widely used to 
' explore the non-perturbative structure of the SDE is to make an ansatz for the full fermion-boson vertex and then 
study the gauge dependence of the physical observables related to the phenomenon of dynamical chiral symmetry 
breaking (DCSB). This method has been quite popular in four dimensional Quantum Electrodynamics (QED). For 
. *-h \ example, the vertex ansatz proposed by Curtis and Pennington, [jl6| , has been extensively used to study the gauge 
dependence of the fermion propagator and the dynamical generation of fermion mass in Quenched QED, e.g., ■ 



Later on, the work of Bashir and Pennington, |21] , |22j| , proposed an improved vertex which achieves complete gauge 
independence of the critical coupling above which mass is dynamically generated. These methods use the cut-off 
regularization to study the gauge dependence of the physical observables. As the cut-off method in general does 
not respect gauge symmetry, a criticism of these works has been raised recently, |23|~|25|]. They suggest dimensional 
regularization scheme to study the chirally asymmetric phase of QED so that the possible gauge dependence coming 
from the inappropriate regulator could be filtered out. 

Three dimensional Quantum Electrodynamics (QED3) provides us with a neat laboratory to study DCSB as 
it is ultraviolet well-behaved and hence the source of gauge non-invariance finds its roots only in the simplifying 
assumptions employed and not in the choice of the regulator. Burden and Roberts, |2(|, studied the gauge dependence 
of the chiral condensate in quenched QED3 and proposed a vertex which appreciably reduces this gauge dependence 
in the range — 1 of the covariant gauge parameter £. Unfortunately, the choice of their vertex does not transform 
correctly under the operation of charge conjugation. Moreover, the selected range of values for£ is very narrow, close 
to the vicinity of the Landau gauge. The studies of QED3 carried out by Bashir et. ai, |^,^8|, reveal that Landau 
gauge may not be a preferred gauge in QED3. In fact, it is seen that the LKF transformation of one loop fermion 
propagator remains insensitive to the inclusion of the LKF transformation of the second loop, provided we are away 
from the Landau gauge. Assuming similar behaviour for the higher loops, we may expect that for larger value of the 
gauge parameter, physical observables would also become insensitive to the value of this parameter. In this paper, 
we undertake the calculation of the Euclidean mass of the fermion (referred to as mass from now onwards) and the 
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condensate for a wide range of values of £. We find that indeed for larger values of the gauge parameter, both the 
quantities seem to become increasingly gauge invariant. 



II. THE FERMION PROPAGATOR 



In quenched QED3, the SDE for the fermion propagator in the Minkowski space can be written as : 

d 3 k 
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where q = k — p, e is the electromagnetic coupling, V(k,p) is the full fermion-photon vertex, S F (p) and A^(q) are 
the bare fermion and photon propagators defined as 



S° F (P) = V 
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and Sp(p) is the full fermion propagator, which we prefer to write in the following most general form : 
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F (p) is referred to as the wavefunction renormalization and A4(p) as the mass function. £ is the usual covariant 
gauge parameter. Eq. (2.1) is a matrix equation. It consists of two independent equation s, w hich can be decoupled 
by taking its trace after multiplying it with 1 and ft, respectively. Making use of Eqs. (2.2,2.3) and replacing the full 
vertex by its bare counterpart, these equations can be written as : 
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Carrying out angular integration after the Wick rotation to the Euclidean space, the above equations acquire the 
form : 
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A trivial solution to Eq. ( |2.7| ) is M.(p) = 0, whic h co rresponds to the usual perturbative solution. We are interested in 
a non-trivial solution by solving Eqs. (2.6) and (2.7) simultaneously. Such a solution for Ai(p) is related to the mass 
m and the chiral condensate < ipip > through the relations m = M.{m) and < iptp >= Ap 2 A4(p) / (2 + , respectively. 
We shall study the gauge dependence of these quantities in the next section for the fixed value of a — 1/4-7T. 



III. EFFECT OF THE WAVEFUNCTION RENORMALIZATION 



In stud ying DCSB, it has been a common practice to make the approximation F(p) = 1 so that we only have to 
solve Eq. (2.7). The justification for this approximation stems from the fact that perturbatively F(p) = 1 + 0(a£/n). 
If a is small and we are sufficiently close to the Landau gauge, one would naturally expect that F(p) « 1. Although, 
it has been quite customary to employ this approximation, there exist several works which include both the equations. 
We study the effects of neglecting the wavefunction renormalization quantitatively. Fig. (Q) depicts the mass function 
Ai(p) for F(p) = 1 in various gauges. As expected, the mass function is roughly a constant for low values of p 
and falls as 1/p 2 for large values of p. The integration region chosen is from 10 -3 to 10 3 and we select 26 points 
per decade. The mass probes low momentum region of this graph, whereas, the condensate is extracted from its 
asymptotic behaviour. Obviously, the mass seems to vary in more or less equally spaced steps with the variation of 



2 



the gauge parameter. In order to obtain a quantitative value of the mass, we select neighbouring points p a and pb 
(p a > pb), such that M{j>a) < Pa and A4(pb) > Pb- We then approximate the mass by the following relation : 

M{pb)~M{p a ) f 

m= (m-p a ) + M(p a ) . (3.1) 

Pb ~Pa 

As for the condensate, the figure does not distinguish between the results for various gauges. Therefore, we have to 
look at the numbers explicitly. Table (1) shows the value of the condensate for £ ranging from 1 — 5. The point 
p = 1000 was chosen to calculate the condensate. This number seems sufficiently large as the 1/p 2 behaviour seems 
to set in much earlier (p ps 300), as noted also in ^6|. In Figs. (|J) and (j|) we display the gauge dependence of the 
chiral condensate and the mass for F(p) = 1 in a wide range of values of the gauge parameter. The condensate varies 
heavily with the change of gauge, roughly twice per unit change in the value of £. Gauge dependence of the mass is 
not too different either. 



£ p ^p 2 M(p 2 ) 


£ p ^p 2 M(p 2 ) 


0.0 1000 2.31109 
642.233 2.31117 
316.228 2.31119 


2.5 1000 11.6992 
642.233 11.7001 
316.228 11.7003 


0.5 1000 3.61103 
642.233 3.61119 
316.228 3.61124 


3.0 1000 14.4432 
642.233 14.4444 
316.228 14.4448 


1.0 1000 5.19982 
642.233 5.20009 
316.228 5.20017 


3.5 1000 17.4761 
642.233 17.4777 
316.228 17.4782 


1.2 1000 5.91622 
642.233 5.91654 
316.228 5.91664 


4.0 1000 20.7977 
642.233 20.7998 
316.228 20.8005 


1.5 1000 7.07745 
642.233 7.07787 
316.228 7.07800 


4.5 1000 24.4081 
642.233 24.4108 
316.228 24.4117 


2.0 1000 9.24390 
642.233 9.24452 
316.228 9.24473 


5.0 1000 28.3073 
642.233 28.3106 
316.228 28.3117 



TABLE 1. The condensate in various gauges for F(p) = 1 



Repeating the exercise by taking both the equations, namely Eqs. ( |2.6| ) and (2.7), into account, we see similar 
qualitative behaviour of the mass function. It is roughly a constant for low values of p and falls as 1 /p 2 for large 
values of p, Fig. (^). The large p behaviour is also evident from the entries in Table (2). As for the wavefunction 
renormalization, it also is constant for small values of p. As p becomes large it goes to 1, Fig. (|). In Table (2) we also 
give a comparison with the work of Burden and Roberts, p6| . As mentioned earlier, they restrict themselves to the 
close vicinity of the Landau gauge, where our results are in excellent agreement. We investigate the gauge dependence 
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of the condensate as well as the mass far beyond the Landau gauge. A graphical description can be found in Figs, (g) 
and (Q). The following points are important to note : 

• The wavcfunction renormalization plays an extremely important role in restoring the gauge invariance of the 
chiral condensate as well as the mass of the fermion. Although the qualitative behaviour of the mass function in 
various regimes of momenta remains largely unchanged, whether or not we employ the approximation F{p) = 1, 
quantitative dependence of the physical observables mentioned above on the covariant gauge parameter £ reduces 
a great deal by including the wavefunction renormalization. 

• As we move away from the Landau gauge towards large positive values of £, the gauge dependence of the 
condensate as well as the mass keeps diminishing, without resorting to any sophisticated ansatze for the fermion- 
boson interaction. 



BHR BR 


£ p ^p2 Mb2) _4_ p 2 M(p2) 

BHR BR 


0.0 1000 2.31109 2.316 
642.233 2.31117 
316.228 2.31119 


2.5 1000 0.98597 
642.233 0.98594 
316.228 0.98571 


0.5 1000 1.77309 1.775 
642.233 1.77313 
316.228 1.77306 


3.0 1000 0.90941 
642.233 0.90938 
316.228 0.90912 


1.0 1000 1.44791 1.447 
642.233 1.44793 
316.228 1.44780 


3.5 1000 0.85201 
642.233 0.85197 
316.228 0.85169 


1.2 1000 1.35288 1.352 
642.233 1.35288 
316.228 1.35274 


4.0 1000 0.80839 
642.233 0.80833 
316.228 0.80803 


1.5 1000 1.23591 
642.233 1.23591 
316.228 1.23574 


4.5 1000 0.77497 
642.233 0.77491 
316.228 0.77458 


2.0 1000 1.09014 
642.233 1.09012 
316.228 1.08992 


5.0 1000 0.74933 
642.233 0.74927 
316.228 0.74891 



TABLE 2. The condensate in various gauges including wavefunction renormalization 



IV. EFFECT OF THE WARD-GREEN- TAKAHASHI IDENTITY 



The bare photon propagator which appears in Eq. fl2.1| ) can be split up in longitudinal and transverse parts as 
follows : 
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where = —g^/q 2 + q^q^/q 4 . Employing this decomposition, we can rewrite Eq. ( [2.l| ) as 



s- F \v) = sl-\v) 



(4.2) 



It is well known that the use of the WGTI, in the equivalent of the last term of Eq. (4.2) in QED4, filters out a spurious 
term which is an artifact of using the gauge dependent cut-off regula tor. Therefore, one is naturally motivated to use 
this decomposition in dimensions other than four. Multiplying Eq. (4^) by 1 and p 1 respectively and Wick-rotating 
to the Euclidean space, we obtain the following equations : 
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As the terms of the type l/(fc 2 — p 2 ) are harder to deal with numerically, we use the approximation F(p) = 1 to 
analyze the effect of the WGTI. Under this simplification, we only have to solve 
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The mass function obtained on solving Eq. (4.7) is depicted in Fig. (|8|), which, along with Table (3), reveals that its 
qualitative behaviour remains unchanged both for small and large values of p. In Figs. (Q) and (|l^) , we compare the 
gauge dependence of the condensate and the mass with and without the usage of the WGTI. As in QED4, we find 
that the gauge dependence of these quantities seems to increase by incorporating the said identity. 
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£ p ^p 2 M(p 2 ) 


£ p ^p 2 M(p 2 ) 


0.0 1000 2.31109 
642.233 2.31117 
316.228 2.31119 


2.5 1000 18.7057 
642.233 18.7085 
316.228 18.7136 


0.5 1000 4.30196 
642.233 4.30221 
316.228 4.30248 


3.0 1000 23.9226 
642.233 23.9268 
316.228 23.9345 


1.0 1000 6.93473 
642.233 6.93529 
316.228 6.93609 


3.5 1000 29.7865 
642.233 29.7924 
316.228 29.8037 


1.2 1000 8.16826 
642.233 8.16900 
316.228 8.17011 


4.0 1000 36.2975 
642.233 36.3056 
316.228 36.3212 


1.5 1000 10.2122 
642.233 10.2133 
316.228 10.2150 


4.5 1000 43.4556 
642.233 43.4663 
316.228 43.4872 


2.0 1000 14.1356 
642.233 14.1375 
316.228 14.1406 


5.0 1000 51.2608 
642.233 51.2746 
316.228 51.3020 



TABLE 3. The condensate in various gauges for F(p) — 1 making use of WGTI 



V. DIMENSIONAL REGULARIZATION METHOD 

In this section we compare our numerical results with those obtained by employing the dimensional regularization 
scheme, [^3|-p5|. For simplicity, we restrict ourselves only to the Landau gauge without incorporating the WGTI. In 
this case, the equation for the mass function acquires the following form in Euclidean space in arbitrary dimensions : 

where a is a dimensionful coupling except in four dimensions. We define d = 4 — 2e and 

a = a d H 2e , (5.2) 

ad being dimensionless. We now use the identity d d k = fc d_1 dk d£ld, where dfld is the d-dimensional solid angle 
defined as dfld — YliZi sin rf_1_i 6i d6i. The angle 9d-\ varies from to 2tt, whereas all other angles vary from to n. 
Choosing 6\ to be the angle between k and p, we can easily carry out the remaining angular integrations to arrive at : 

(47r)Vr(^i) io fc 2 + A4 2 (p) J q 2 
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Using the standard formula, p9| , 
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This equation was studied in detail in [g3| in four dimensions, taking e to be a small positive number. The factor 
(fc 2 U e in the numerator regulates the otherwise divergent behaviour of the integrand for large momenta. As noted 
in p3], the hypergeometric function does not play any role in regularization and hence can simply be replaced by 
F(l,0;2,z) = 1. In case of three dimensions, the hypergeometric function develops a pole for k 2 — p 2 , as is obvious 
from the following identity : 



F(l,e;2-e;l) = 



1 



1 - 2e 



(5.5) 



As was pointed out earlier, such terms are hard to deal with numerically. Due to increasing computational time and 
memory, we go only up to e = 0.48, starting from e = 0.4. To obtain satisfying results, we need to use increasingly 
more points per decade as we approach closer to e = 0.5. For instance, we use 100 points per decade for e = 0.46. 
The problems of ever increasing computational time and memory limited us to use 140 points per decade for the case 
of e = 0.48. Despite this large number, we belive that the corresponding result falls short of the desired accuracy. As 
a result, there is a slight rise at the end of the fiat regi on o f the mass function, and the final descent begins rather 
late, Fig. ([n]). The problematic pole for e = 0.5 in Eq. ( |5.5[ ) corresponds to the relatively well-controlled singularity 
in the following expression 



F 



1 3 2 

'2'2' Z 



1 . 1 + 2! 

= — In 

2z 1-z 



(5.6) 



for z — ► 1. A comparison between the mass function obtained from techniques based upon the dimensional regular- 
ization scheme and the one computed in Section III is also depicted in Fig. ([ll]). Taking the numerical limitation for 
e = 0.48 into account, we note that as e approaches the value of 0.5, we get closer and closer to the result obtained 
in Section III, where we work in 3-dimensions to start with. 



VI. CONCLUSIONS 



The sources of gauge non-invariance in the study of SDE in QED can arise from (i) the use of the approximation 
F(p) = 1, (ii) the violation of WGTI and the LKF transformations, (iii) the flawed ansatz for the vertex and (iv) 
the inadequate choice of the regulator. In case of 3 dimensions, QED becomes neater because of the absence of 
ultraviolet divergences. The regulator dependence does not pose a threat any longer. In this paper, we have studied 
the quantitative effect of the approximation F(p) = 1 and the violation of the WGTI on the gauge dependence of 
the mass and the condensate in QED3 for a wide range of values of the covariant gauge parameter £. Lifting the 
approximation F(p) — 1 seems to reinstate gauge invariance to an impressive extent. Partial use of the WGTI seems to 
increase the gauge dependence, an effect noted also in 4 dimensions in |23| . As an advantage of using a broad range of 
values of £ , we find that when they are large, the gauge dependence of the said quantities diminishes without resorting 
to improved ansatz for the vertex. Currently, work is underway to explore the role played by various elaborate choices 
of the vertex in this context. 
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Mass Function in Various Gauges 

(Bare vertex with F=l and without the use of WGTI) 
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FIG. 1. The mass function M(p) in the approximation F(p) = 1. 



Condensate in Various Gauges 

(Bare vertex with F=l and without the use of WGTI) 



FIG. 2. The condensate < ipip > in the approximation F(p) = 1. 
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Mass in Various Gauges 

(Bare vertex without the use ofWGTI with F=l) 



FIG. 3. The mass in the approximation F(p) = 1. 
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Mass Function in Various Gauges 

(Bare vertex without the use of WGTI) 
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FIG. 4. The mass function including the equation for the wavefunction renormalization. 
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FIG. 5. The wavefunction renormalization. 
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Condensate in Various Gauges 

(Bare vertex without the use of WGTI) 



F=l 
including F 



FIG. 6. The condensate: a comparison between the cases with and without the use of F(p) = 1. 
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FIG. 7. The mass: a comparison between the cases with and without the use of F(p) = 1. 



Mass Function in Various Gauges 

(Bare vertex with F=l and with the use of WGTI) 
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FIG. 8. The mass function including the effect of the WGTI for F(p) = 1. 
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Condensate in Various Gauges 

(Bare vertex with and without the use of WGTI for F=l) 
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FIG. 9. Effect of the WGTI on the gauge dependence of the condensate for F(p) = 1. 



Mass in Various Gauges noWGTi * 

(Bare vertex with and without the use of WGTI for F=l) 



FIG. 10. Effect of the WGTI on the gauge dependence of the mass for F(p) = 1. 
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Dimensionally Regularized Mass Function 42 
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FIG. 11. The mass function for various values of e. Note that e = 0.5 would correspond to 3 dimensions. The result for 
e = 0.48 could not be achieved up to the desired accuracy because of the ever increasing computational time and memory 
involved. 
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